A PROOF OF THE POSITIVE DENSITY CONJECTURE FOR INTEGER APOLLONIAN 

CIRCLE PACKINGS 



JEAN BOURGAIN AND ELENA FUCHS 

Abstract. A bounded Apollonian circle packing (ACP) is an ancient Greek construction which is 
made by repeatedly inscribing circles into the triangular interstices in a Descartes configuration of four 
mutually tangent circles. Remarkably, if the original four circles have integer curvature, all of the circles 
in the packing will have integer curvature as well. In this paper, we compute a lower bound for the 
number k(P,X) of integers less than X occurring as curvatures in a bounded integer ACP P, and prove a 
conjecture of Graham, Lagarias, Mallows, Wilkes, and Yan that the ratio k{P,X)/X is greater than for 
X tending to infinity. 



1. Introduction 

In the first picture in Figure[T|there are three mutually tangent circles packed in a large circle on the 
outside, with four curvilinear triangles inbetween. By an old theorem (circa 200 BC) of Apollonius of 
Perga, there are precisely two circles tangent to all of the circles in a triple of mutually tangent circles. 
One can therefore inscribe a unique circle into each curvilinear triangle as in the second picture in 
Figure [T] Since this new picture has many new curvilinear triangles, we can continue packing circles 
in this way - this process continues indefinitely, and we thus get an infinite packing of circles known 
as the Apollonian circle packing (ACP). Since the radii of most of the circles in an ACP are extremely 
small, it is convenient to measure the circles in the packing via their curvatures, or the reciprocals of 
the radii, instead. 




Figure 1 . Apollonian Circle Packings 
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A remarkable feature of these packings is that, given a packing in which the initial four mutually 
tangent circles have integer curvature, all of the circles in the packing will have integer curvature as 
well. Such a packing is called an integer ACP and is illustrated in Figure|2]with the packing generated 
by starting with circles of curvatures 1,2,2, and 3. 

There are various problems associated to the diophantine properties of integer ACP's which are 
addressed in IGLMW Y I by the five authors Graham, Lagarias, Mallows, Wilks, and Yan. They make 
considerable progress in treating the problem, and ask several fundamental questions many of which 
are now solved and discussed further in ED, E, 113, llKOl , and llS2l . 

In most of these papers, ACP's are studied using a convenient representation of the curvatures 
appearing in an ACP as maximum-norms of vectors in an orbit of a specific subgroup A of the orthog- 
onal group 0(3, 1). We introduce this group in Section [l~T| and will use it throughout. In regards to 
counting the number of integers represented in a given ACP, Graham et. al. exploit the existence of 
unipotent elements of A in OGLMWYII to establish the lower bound below for the number k{P,X) of 
distinct curvatures less than X of circles in an integer packing P: 

(1.1) k{p,x):$>Vx 

where the notation 

y >/3 z or y <|3 z 

in this paper is taken to mean that there exists a constant c > depending only on j3 such that 

y>cz or, respectively y < cz- 

Graham et. al. suggest in IGLMWYI that the lower bound in can be improved. In fact, they 
conjecture that the integers represented as curvatures in a given ACP actually make up a positive 
fraction of the positive integers N. 

It is important to note that this question is different from one recently addressed in IKOII by Kon- 
torovich and Oh about the number Np{X) of circles in a given packing P of curvature less than X. 
This involves counting curvatures appearing in a packing with multiplicity, rather than counting every 
integer which comes up exactly once as we do in this paper In fact, the results in IIKOI suggest that 
the integers occurring as curvatures in a given ACP arise with significant multiplicity. Specifically, 
Kontorovich and Oh find that Np{X) is asymptotic to c-x^, where 5 = 1.3056... is the Hausdorff 
dimension of the limit set of the packing. Kontorovich and Oh's techniques, however, do not extend in 
any obvious way to proving that the integers represented by curvatures in an ACP make up a positive 
fraction in N. 
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In BSII Samak uses the existence of arithmetic Fuchsian subgroups of A to get a bound of 

X 



(1.2) 



k{P,X) > 



towards Graham.et.al.'s positive density conjecture. This method, which we summarize in Section|2] 
was further improved to yield a bound of 

X 

k(RX) > 

where e = 0. 150 ... by the second author in a preprint IIFll . 




Figure 2. Apollonian Circle Packing (-1,2,2,3) 



In this paper, we refine this Fuchsian subgroup method in a number of ways and settle the positive 
density question of Graham et.al. in the following theorem: 

Theorem 1.1. For an integer Apollonian circle packing P, let k{P,X) denote the number of distinct 
integers up to X occurring as curvatures in the packing. Then for X large we have 

k{P,X)^X 

where the implied constant depends on the packing P. 

We treat this question by counting curvatures in different "subpackings" of an ACP. Namely, we 
fix a circle Co of curvature oq and investigate which integers occur as curvatures of circles tangent 
to C. This gives the preliminary lower bound in \l-2\ which was first proven by Sarnak in IS II . The 
essential observation which leads to this lower bound is that the set of integers appearing as curvatures 
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of circles tangent to Co contain the integers represented by an inhomogeneous binary quadratic form 

/flo(^'3')-flo 

of discriminant —A-Oq. Our approach in Section [s] is to repeat this method for a subset of the circles 
which we find are tangent to Co in this way. For every circle C of curvature a tangent to Cq, we can 
produce a shifted binary quadratic form 

fa{x,y)-a 

where /„ has discriminant — 4fl^ and consider the integers represented by /„. We consider a in a 
suitably reduced subset of [(logX)^, (logX)^] and count the integers represented by /„ — a for a in 
this subset. It is important to note that the integers represented by /„ and /„/ for d are a subset 
of integers which can be written as a sum of two squares since both forms have discriminant of the 
form —5^. In fact, /„ and /„/ represent practically the same integers (see IFlll for a more detailed 
discussion). It is rather the shift of each form /« by a that makes the integers found in this way vary 
significantly. Our final step is to give an upper bound on the number of integers in the intersection 

{m represented by — a} n \m represented by /„/ — a'} 

In obtaining this upper bound, we count integers with multiplicity, which is a sacrifice we can afford 
to make for our purposes. This method leads to a proof of the conjecture of Graham et.al. that the 
integers appearing as cui-vatures in a given integer ACP make up a positive fraction of all integers. 

Acknowledgements: We thank Peter Sarnak and Jeffrey Lagarias for introducing this problem to 
us, and for many insightful comments and conversations. We thank Alex Kontorovich for sharing his 
program for drawing Apollonian circle packings. 

1.1. The Apollonian group. In 1643, Descartes discovered that the curvatures of any four externally 
cotangent circles of curvatures x\, xi, x^,, and x\ satisfy the equation 

(1.3) g(xi,X2,X3,X4) = 2{x\ +X2 +X3 +X4) — (xi +X2 +X3 +X4)^. 

A proof of this can be found in OCxI . This equation also holds for four mutually tangent circles 
when three of the circles are internally tangent (or inscribed) to the fourth circle as in the case of 
ACP's. However, because the outside circle is internally tangent to the other three, we assign to it 
a negative curvature in order for the equation to hold. We note that Descartes' theorem allows us to 
solve a quadratic equation for the fourth curvature once three are known (this is one way to verify 
ApoUonius's theorem about circles tangent to a triple of mutually tangent circles). 

We thus assign to every set of 4 mutually tangent circles in a packing P a vector v G of the 
circles' curvatures, and use Descartes' equation to express any ACP as the orbit of a subgroup A of 
the orthogonal group Og(Z) acting on v. In IGLMWY I. Graham et. al. describe an algorithm for 
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generating a root quadruple v, which corresponds to the four largest circles in the packing, and is 
uniquely determined for any given ACP. From this point on, we denote by v the root quadruple of the 
packing in question, and assume that the packing P is primitive - the gcd of the coordinates of any 
w G A\ in the orbit is 1. The group A, called the Apollonian group, was introduced by Hirst in 1967 
(see IHl). It is a group on the four generators 
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derived by fixing all but one of the x, in (1 .3 1 and solving Q{xi) — for the fourth. Note that each 5, is 
of order 2 and determinant —1, and that left multiplication of 5/ by a vector w = (wi ,W2, W3, W4)^ g 
fixes three of the coordinates of w. We denote — Ay the orbit of A acting on v, and note that our 
results hold for an arbitrary v, as they are based on properties of the group A, and do not depend on the 
specific orbit in question. Since the quadratic form in ( 1.3 1 is of signature (3,1) over M, we have that 



A is a subgroup of 0(3, 1) and can be thought of as a subgroup of the group of motions of hyperbolic 
3-space H^. In this way A is a discrete algebraic group acting on where the complement of three 
mutually tangent hemispheres inside an infinite cylinder is the fundamental domain of the action. This 
fundamental domain has infinite volume, which makes counting integers in the group's orbit quite 
difficult. We note, however, that A contains Fuchsian triangle subgroups generated by any three of the 
Si above, which are lattices in the corresponding 0(2, l)'s. We use this fact extensively throughout 
this paper. 

To this end, denote by A, the subgroup of A generated by three of the four generators as follows: 

A,:=({5i,52,53,54}-{5i}). 

This group is the Schottky group generated by reflections in the three circles intersecting the ith circle 
in the root quadruple and perpendicular to the initial circles in the packing; in particular, the ;th circle 
is fixed under this action. The fundamental domain of A,- is then a triangle bounded by the three circles, 
and has hyperbolic area %. 
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2. A PRELIMINARY LOWER BOUND 

In this section, we follow MSII in order to count integer points in an orbit of a subgroup A, of the 
Apollonian group as described in Section [T7T| This produces a preliminary lower bound on the number 
k{P,X) of integers less than X occurring as curvatures in an Apollonian packing P. 

Proposition 2.1. For an integer Apollonian circle packing P, let k{P,X) denote the number of distinct 
integers less than X occurring as curvatures in the packing. Then we have 

(C(P,X)»-^. 

VlogX 

Proof. We fix a circle Co of non-zero curvature aq in the packing P, and count the integers which occur 
as curvatures of circles tangent to Co in P. This is identical to considering the orbit of Ai acting on a 
quadruple v of mutually tangent circles {ao,b,c,d), since A i fixes the first coordinate of v and its orbit 
represents all of the circles tangent to C„. Note that A generates all possible Descartes configurations 
in the packing P, and there can only be finitely many circles of curvature aq in the packing since 
the total area of all the inscribed circles is bounded by the area of the outside circle. Therefore it is 
reasonable to count the circles represented in the orbit of Ai, since they make up a positive fraction of 
all of the circles in P tangent to a circle of curvature aQ. 

In this orbit, we have that the first coordinate ao is fixed, and the other coordinates of points in the 
orbit of Ai vary to satisfy 

Q{ao,X2,X3,X4) = 2{al +X2 +x\ +X4) — (ao +X2 +X3 +^4)^ ~ 0, 

where Q is the Descartes form in ( |1.3| l. A change of variables y = (3'2,}'3,}'4) = (.*^2j.*^3 7'^4) + (flo,flOjflo) 
allows us to rewrite the equation above as 

(2.1) g{y)+Aal = 0, 

where g{y) —y\+y\^y\~ lyiy^ — '^yiyA — '^y^yA is the resulting ternary quadratic form. We can thus 
conjugate the action of Aj on (flo,X2,X3,X4) to an action independent of ao which preserves the form 
g. This is the action of a group F on y, generated by 

1 

2 -1 

Moreover, the action of F on 

v' = (/j + aojC + aoj'^ + oo) 

is related to the action of Ai on v by 

Aiv= (ao,r[v'-(flo,flo,flo)]), 
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SO we count the same number of curvatures occurring in the packing before and after this change of 
variables. We change variables once again by letting 

y2=A,y3^A + C~2B, y^ = C. 

We note that (y2,3'3,3'4) G 2^ implies that A, B, and C are integers, and the primitivity of the packing 
is preserved as well - the gcd of A,B, and C is 1. With this change of variables, Y is conjugated to an 
action of a group Y' on (A,B,C) which is generated by 
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Under this change of variables, the expression in ( |2.1[ i becomes 



(2.2) 



4(^2 -AC) 



Letting A(A,B,C) denote the discriminant of the binary quadratic form Ax^ + 2Bxy + Cy^, (2.2 1 is 
simply 

A{A,B,C)^al, 

and thus Y' is a subgroup of Oa(Z), the orthogonal group preserving A. Let Y denote the intersection 
r' n SOa(Z). The spin double cover of SOa is well known (see BEGMI ) to be PGL2, and is obtained 
via the homomorphism 



p :PGL2(Z) 



SOa(Z) 



(2.3) 



a j3 
7 5 



aS-Py 



a" 



lay 



aj5 ad + jSy y5 
V 2^5 



f ) 



written here over Z as this is the situation we work with. It is natural to ask for the preimage of f 
under p which we determine in the following lemma. 

Lemma 2.2. Let Y and p he as before. Let A(2) he the congruence 2-subgroup ofPSL2{'Z). Then the 
preimage ofY in f GL2(Z) under p is A(2). 



and 



Proof. We can extract from the generators of Y' as well as the formula in ( 2.3 1 that the preimage of f 
under p contains 

: : ::-'). 

and so f contains the congruence subgroup A(2) of SL2(Z). Recall that the area of Ai\IHI^ is 
and note that SOa(Z) HT' contains exactly those elements of Y' which have even word length when 
written via the generators of F', making up half of the whole group. Therefore the area off\H2 is27i;. 



1 

-2 1 
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which is equal to the area of A(2)\H2, and hence the preimage of f in PGL2(Z) is precisely A(2) as 
desired. □ 

Recall that we would like to count the integer values of y2,y3, and 3^4 - in terms of the action of F, 
we are interested in the set of values A,C, and A + C — 2B above. Lemma [2!2] implies that these values 
contain integers represented by the binary quadratic form 

(2.4) /„„(C,v)=AoC'+2BoCv + Cov2, 

where {^,v) — 1, and the coefficients are derived from the change of variables above: 

b + d-c 

(2.5) Ao=b + aQ, Co = d + ao, Bo = . 

We note that the discriminant of this form is not a square, since 

(2.6) D(/„o) = (2Bo)^-4AoCo = -4flg 

and flo ^ 0. Since the vectors in the orbit of Ai are of the form {ao,A — ao,A + C — 2B — ao,C — ao), 
they correspond to the integer values of 

(2.7) fag{x,y)-ao, 
where is as before. Therefore 

(2.8) ^iP,X) > #{m e Z |m > 0, fag{x,y)-ao = m for somex,y € Z, {x,y) = 1} 

and we need only to count the integers represented by fag in order to get a bound on the number of 
curvatures in P. This is done in | J |, from which we have the following: 

Lemma 2.3. (James): Let f be a positive definite binary quadratic form over Z of discriminant —D, 
where D is a positive integer. Denote by Bd{X) the number of integers less than X represented by f. 
Then 

cX ( X 
Bd(X) = +01 

where 



?=3 (4) 



Lemma [273] paired with (2.8 1 implies that 



KiP.X) > , 

as desired. □ 
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3. Proof of Theorem I1.1I 



In this section we sharpen the bound in Proposition 2.1 in order to answer the question posed by 
Graham et.al. in P GLMWYl and prove that the integers appearing as curvatures in any integer ACP 
make up a positive fraction of all positive integers. Our computation in Section |2]reflects only those 
circles which are tangent to a fixed circle in P. It is thus natural to count some of the omitted curvatures 
here. Specifically, we repeat the method from Section |2] several times, fixing a different circle C each 
time and counting the integers occuring as curvatures of circles tangent to C. 

Recall that to prove Proposition |2. 1 | we fixed a circle of curvature oq, and associated curvatures of 
circles tangent to it with the set of integers (without multiplicity) represented by fag{x,y) — ao- We 
denote the set of these integers that are less than X by jz^q: 

— {a EN\a < X,faQ{x,y) — ao = a for some integers x,y > 0} 

For every a € we use the method in Section|2]to produce another shifted binary quadratic form 

f„{x,y)-a 

of discriminant D = — 4a^. As in Section [5] we wish to count the integers represented by these new 
forms. For each a e x/q, let Sa denote the set of integers less than X represented by /„ — a: 

Sa = {n E N\n < X ,n — fa{x,y) — a for some relatively prime integers x,y >Q} 

Note that the sets S^, depend only on oq, the curvature of Cq. One important consideration in counting 
the integers represented by the forms /„ is that their discriminants can be very large with respect to 
X, and thus many of the represented integers may be > X. In particular, the count in Lemma 2.3 is 
not uniform in D so we use more recent results of Blomer and Granville in liBG | which specify how 
the number of integers less than X represented by a binary quadratic form depends on the size of the 
discriminant of the 

With this notation, the bounds in HBGI yield a lower bound on Y,a \Sa \ for the a's we consider. We 
also compute an upper bound on a'\^a^ S^,/ \ for a ^ a' so that 

^|5„|-£|5„n5„,| 

a a.a' 

gives a lower bound for k{P,X). A crucial ingredient to computing this and proving Theorem ] 1.1 [ (that 
the integers appearing as curvatures in a given ACP make up a positive fraction of N) is the balance 
between these lower and upper bounds - for example, the more sets Sa we choose to include in our 
count, the bigger the lower bound on Y.a \Sa\- However, choosing too many such sets will also increase 



2 

The results of Blomer and Granville concern quadratic forms of square-free determinant, but the authors note in section 9.3 of 
IBGI that the same can be done for binary quadratic forms of non-fundamental determinant as well. 
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the upper bound on the second sum Y.a.a' \Sa n5„'|. In fact, it is possible to choose so many sets Sa 
that the upper bound on the intersections outweighs the lower bound on the sizes of Sa- In SectionjTT] 
we specify how we choose the a's used in our computation, and compute the first sum, ^„ \Sci\- In 
we compute an upper bound on ^„ „/ \Sa H 5„/ 1 for a 7^ a' to prove Theorem ] 



Section 



3.2 



3.1. Integers represented by multiple binary quadratic forms. In this section, we evaluate the 
sum Y,a \Sa\, choosing a's in a subset of m/q in order to ensure that we obtain a positive fraction of X 
in our final count. Specifically, we consider a e M) such that 

{\ogXf <a< (logX)^ 

This interval is chosen to give us the desired lower bounds in conjunction with results in fBGl - this 



will become clear in the computations preceding ( 3.9 1. We would like to further reduce the set of a's 
we consider so that the bounds on the size of the intersections of sets Sa are not too large. To do this, 
we first partition the interval [(logX)^, (logX)^] into dyadic ranges [2'^, 2'''+'] and select a's within 
these ranges. 

Namely, we consider ^0 n [2^,2*^+'] where {logXf < 2'', 2''+^ < {logXf. The size of this set 
depends only on ao, the curvature of the original circle we fixed. By Lemma [23] we have 



(3.1) 



k ^k+l- 



> 



into intervals 



where the implied constant depends on ao. We partition each dyadic interval [2 ,2 
[2*^ + n • 77-^,2* + (« + 1) • 77-^] of length T] -^, where < 77 < 1 is a fixed parameter whose im- 

We note that the average over < n < \/krj ^ ' of 



3.3 



portance will become apparent in Proposition 
cardinalities of the corresponding subsets of jz/q is 

(3.2) E, 



{,n[2^- + «-Tj^,2'^ + (« + i).T7^: 



by ( 3. 1 1. Thus for every value of k there exists an < « < Vkrj ^ for which the intersection in ( 3.2 1 



contains ^ J] ^ integers. For simplicity of notation, we assume without loss of generalit}] |that n — 0, 
and define to be 

2' 



(3.3) 
where we have 
(3.4) 



=£/on[2^2'^- 



2' 



One can in fact extend Lemma 



2.3 



14.4 of IfTI . However, it is not necessary liere. 



to show that this holds for every n. Friedlander and Iwaniec do this for ao = 1 in Theorem 
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up to a constant which depends only on ciq. Denote the union of these subsets by £/ : 

(3.5) ^ = 1J^W 

The resuhs in LEG J imply the following lemma regarding the integers represented by quadratic forms 
associated with a£ £/. 

Lemma 3.1. Let and fa be as before. Then we have 

To prove Lemma |3.1[ we recall the notation and relevant theorem from OBGl . Let / be a binary 
quadratic form of discriminant D, and let r/(n) be the number of representations of n by /: 

(3.6) r/(n) =#{(oti,OT2) G - {0} | gcd{m\,m2) = \ ,f{m\,m2) = n} 

Let U'^{X) — Ln<x the number of integers less than X represented by /, counting without 
multiplicity. In OBGl . Blomer and Granville compute bounds for U^{X) for D in three ranges between 
and X. These ranges are defined in terms of the class number h of the binary quadratic form / and 
by g, the number of genera. Letting £ = = L{\ ,X-d){^{D) / D), they create a parameter 

log(/z/g) 
(log2)(log(^_DlogX)) 

where h/g = Their bounds for U^f{X) are then uniform in D for each range below (see 

Lemma 13. 2b: 



• < K- < 1/2 

• 1/2 < K" < 1 

— ^ log log £> 



In the first and last range, they are able to compute both an upper and lower bound on U . However 

the smallest range is essentially James' result in Lemma [23] and is used to show that 

X 



they prove only an upper bound for U^{X) in the case that D is in the middle range, which is not 
suitable for our purposes. The lower bound for U%X) for a form / of discriminant —D where D is in 



k{P,X) > 



{\ogxY 

in OFIL Our results and the statement in Lemma [TT] depend on Blomer and Granville's lower bound 
for Uj{X) where / is of dicriminant —D and D is in the third range above. Specifically, we use 
Theorem 2 from IBGI . which is summarized in the lemma below. 

Lemma 3.2. ( Blomer, Granville ): Let f be a binary quadratic form of discriminant —D, and let t/j? {X) 
be as before. Let be the group of genera of binary quadratic forms of discriminant — D. Denote by s 
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the smallest positive integer that is represented by f, and by u the smallest positive integer represented 
by some form in the coset fiS . Then 



(3.7) U%X)^%-[\-^\-^^E^(X,D) 



where 

fx , , (X\QgX X 
(3.8) £o(X,0)<-^- + t(D)-(^— ^ + 

where t{D) is the number of prime divisors ofD, and the implied constant does not depend on D. 



With this in mind we are ready to prove Lemma 3. 1 
Proof of Lemma \3.1\ 

We use Lemma |3.2| to count the integers less than X represented by forms of discriminant D 
where D is a power of \ogX in our case. Recall that 

fa{x,y)^ax^ + 2^xy + yy^ 

is of discriminant D = —Aa^. In particular, since [XogX)^ < a < (logX)^, we have that 

{logX)" <D<{logX)^ 

and the number of prime divisors of D is 

T(D)<loglogX, 

and so 

^°(^'^)«(1^ + (1°^^)-^ 



Thus we have that the error Eq{X,D) <C 3^,^ for any £ > 0, and thus Lemma 3.2 implies 



£.3/4 

(3.9) U%X) > ^ 

\/D 



where the implied constant does not depend on D. Since D — — 4a , it follows from (3.9 1 that the 
number of distinct values less than X represented by fa is S> f and we have 

EN » 

» X E \ 

2*'<(logX)3 
2*>(logX)2 

(3.10) > 77 

as desired. □ 
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The sum in Lemma [3T| is the lower bound on the number of integers we count by considering the 
quadratic forms associated with a G £/. In order to prove Theorem |1.1[ we obtain an upper bound on 
the number of integers we have counted twice in this way in the next section. 

3.2. Integers in the intersections. To prove Theorem |1.1| we would like to show 



(3.11) 



since this union is a subset of all curvatures less than X in the packing P. By Lemma 3.1 we may 
estimate the size of this union as follows: 



(3.12) 



We need only to determine an upper bound for the last sum above. We do this by counting points 
{x,y,x'y') in a box on the quadric 

fa{x,y) - fa'{x ,y') ^ a' - a 

for each a ^ a' £ . The region in which we count these points is induced by the condition that 
/a(jc,y) < X. Namely, rewriting the binary form /„ as 

(ax + j3y)2+4flV 



(3.13) 
we can define a region 



fa{x,y) 



a 



(3.14) 



Ba = {{x.y) e M.^ s.t. \ax + l5y\ < V|a| and \y\ < -^^^} 



a 



so that fa{x,y) < 1 for {x,y) e B„, and fa{x,y) < X for every {x,y) G \/XBa as desired. Therefore, 
the region in M^* over which we consider the forms /„ — /„/ will be 

;^a.a'^{^Ba-^VXBj)r\t^ 

With this notation, we are ready to prove the following proposition. 

Proposition 3.3. Let si , Sa, Sj, rj, andX be as before. Then there exists c > depending only on oq 
such that 



(3.15) 



^ \SanSj<cri^X 



Note that, since we chose < 7] < 1, we have 7] < Tj, and so this upper bound on the size of the 



intersection of the sets 5„ — a is small compared to the count in Lemma 3.1 
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Proof. We note that the expression inside the sum has an upper bound 

(3.16) \SaC^Sa>\ 

< \{{x,y,x' ,y') e \fa{x,y)-fa'{x',y') =0-0'}! 



Although bounding ( 3. 15 1 in this way involves counting the integers in Sa H 8^,' with multiplicity, our 
analysis shows that this sacrifice is in fact not too expensive to our final count. We thus consider the 
quaternary quadratic fom 

F{x,y,x',y') = fa{x,y) - fa'{x' ,y') 

with discriminant A = ()3^ — a'y){^''^ — o:'/) = I6a^{a')^. To obtain an upper bound on the number 
of points in X G ^a.a' for which F{x) =a' — a, one can use the well developed circle method following 
Kloosterman in IJK) and Esterman in [El or modular forms (see MDRSI ). Both methods would yield 
what we want - the latter would give the best results but is not as flexible as the former for our purposes 
since we wish to vary the parameters a and a' which is more straightforward in the circle method. 
Heath-Brown's Theorem 4 in |HB | and Niedermowwe's Theorem 5.6 in |N | determine representation 
numbers of a fixed indefinite quadratic fornj^ Since our a, a' are all a small power of \ogX, the proofs 
of these theorems can be manipulated slightly to yield the following lemma regarding representation 
numbers of all the indefinite quaternary quadratic forms we consider: 

Lemma 3.4. Let F be as before, and let (logX)^ ^a^a' < iXogX)^. Let Xa,a' denote the characteristic 
function on the region i^a,a'' cind let 

F{x)=a—a' 

Let A be as above. Then we have 

Ry ,(a — a') 

f X ■ A^™ \ 

(3.17) = |/^^^,(«_,')|.|e(a-a')|+0(^^j^j 

where the first factor is the singular integral 



(3.18) /; 



7 / 



XaM'{x)e{z{F{x)-a + a'))dx 



dz 



and the second factor is the singular series 
(3.19) 6(fl-a')=n^/ 



p 

p 



^Note that in (HB] one considers representations of an integer m by F where m is asymptotic to the scaling factor P of the 
unsealed domain 6 (in our case B = BaX B^i and P = \/X), while in |N| m is any nonzero integer. 
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where 
(3.20) 



Op = lim p-^^ ■ #{x e {^/pif s.t. F{x) = a-a' (p*)} 



and e(z) = e 



— ^2niz 



In the error term in (3.17 1, A is an arbitrary large fixed constant. With more effort we can in fact 
get a power saving here by using modular forms - while this would yield the best result, the methods 
in BHBII and IjN] suffice. In particular, the argument in IN] lends itself well to our consideration of 
the quadratic form F, which has a discriminant of size (logX)''. The error term in Niedermowwe's 
Theorem 5.6 consists of a power saving in X, and a careful examination of the proof shows that the 
dependency on the discriminant of the form is absorbed into the error term since it is only logarithmi- 



cally large - this is reflected in (3.17 1 via a power of the discriminant A of F. It is similarly important 



here that the distorsion of Ba and B„/ with respect to the standard cube discussed in MHBI and ||Nl is 
logarithmic in X. 

To prove Proposition |3.3| it remains to evaluate the singular integral and singular series in p.l8| l 
and (3. 19 1. For a set P C M^, let V(f ) denote the measure of P. From the definition of /„ and B^, in 



( |TT3| l and (|TT4]), we have 



(3.21) 



lim -■Y({{x,y,x',y') E VXB„ x VXB,, \ \fa{x,y) ^ f„,{x' ,y') - a + a'\ < e} 



lim - 

e^O £ 

X 



X ^/\¥\x 



a' 



To evaluate the singular series (5 (a — a') we prove the following lemma. 



Lemma 3.5. Let & (a~ a') be the singular series defined in (3.19 L We have 



6(fl-fl')« n (i+-)-2™«"-"'» 



p\aa^{a—a') 



Proof. We compute an upper bound for the expression in the limit in (3.20 1 by letting k ~ 1 since the 
expression in the limit decreases with k. Note that if p\{a,a'), we have F is not degenerate modulo p 



by the primitivity of the packing and the definition of the coefficients of /„ in (2.5 i. Therefore 

dp < •#{x e (^/pz)4 s.t. F{x) = (p)}, 

and over Fp, the number of nontrivial representations of by F is bounded above by 2p^ (see (Cl, 
for example), so Gp is bounded above by 2 in this case. In the other cases, we use exponential sum 
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estimates. Taking A: = 1 as before, we have 



P r=0 x,y x'.y 

where ep{z) — exp(^^). There are several cases to consider: 



Case 1: p does not divide aa'{a — a'): 
If we diagonalize /„ and /„/, we obtain 

(^p=p + -L I - 

P r= 1 V P 



\ 2 / \ 2 

ar\ I (x'r\ 2 



p-ep{r{a-a'))^ p +o{p) 



since {a — a',p) = 1. 

Case 2: p\a — a' and does not divide aa': 
In this case we have Gp < p^ +o{p^). 



Case 3: p\a and p Jo': 

Diagonalizing f^i, we obtain 



ar 



. ■ L ( ~~ I ^'v^'/''^P <f^+o(/'^) 

P r=l \P J 

From these bounds and Lemma IsTtI we obtain the desired result in Lemma [331 



□ 



Combining our computation of the singular integral in ( |3.21[ i and the bound on the singular series 
in Lemma [33] the result of Niedermowwe in Lemma [J4| yields 



(3.22) 



|5„n5„,|«^- n (i + -)-2»«"'"') 

aa , , , V 



p\aa'[a—a') 

p\{a,a') 



where (X)(n) is the number of distinct prime factors of n. Thus to evaluate the last sum in (3.12 1, we 
count the number of a G ^ in progressions a = r mod q. To this end, we recall Theorem 14.5 from 
IIfH of Friedlander and Iwaniec regarding sums of squares in progressions in the following lemma 



'''Note that the set of integers in the interval [2*, 2* + 7] ^] which can be written as sums of two squares contains the a G ^e/'*' 

in progressions a = r[q), since jz/'*' is a set of integers represented by a binary quadratic form of discriminant —5^. This count 
is therefore an upper bound on what we want. 
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Lemma 3.6. (Friedlander, Iwaniec): Let b[n) be a characteristic function defined as 

bin): 

and let 



1 ifn — s'^+t^forsomeSjt^'E 
otherwise 



B{x,q,a):^ Y, b{n) 

n<x 

n=a(q) 

For 2 <q <x, (a.q = I), and a = I mod {4-,q) we have 



B{x,q,a) = — ■ 



q ^/\ogx 



1+0 



1 \ n 

logx 



where the implied constant is absolute and Cq <^ log log ^ is a positive constant. 



We note that the statement in Lemma 3.6 is much stronger than what we need - we require only 
an upper bound on B{x,q,a), which could be proven using an upper bound sieve. Since our set s/ 
is obtained via the fixed quadratic form of discriminant — 4aQ from Section |2] such an upper bound 
implies the following in our case. 

Lemma 3.7. Let , X, and rj be as before. Then we have 

l^^loglog,.^ 

a=r{q) 

where \ <q< logX is a square-free integer. 



Proof With the definition of si^^"^ in ( 3.3 1, we may bound above the sum in Lemma 
k for which (logX)^ < 2*, 2*^+1 < (logX)^: 



3.7 



as a sum over 



(3.23) 



a6[2*,2«:+J} 



By Lemma 3.6 the inner sum is bounded above (up to a constant) by 



■n 



c^2^ 
q k 



l+O 



log log ^ 



Since Cq <C loglog^, substituting this into (3.23 1 we have 
as desired. 



log log?. 



< 7? 



log log g 

q 



□ 
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With this in mind, we may evaluate the sum in p.22| i as follows. 



(3.24) 



(3.25) 



£ \Sans,>\ « £ —T 



<«((«,«')) 



aa 



n 1 



p\aa'{a—a') 

p\(a.,a') 



10,11,11,12 



2(0(12) 



E - 

aa' 

loi'iW 



where qQ,qi,q\,q2 are square-free and relatively prime. We may restrict to primes p < (logX)™ in 
the product in (3.24 1, we may restrict in (3.25 1 the summation to qo,qi,q'i,q2 < (log^)™. We bound 
the sum 

L - 

aa' 

101\ \'i 

loi'iW 

121"-"' 



using Lemma 3.7 First fix a and sum over a' subject to the restrictions qoq'i\a' and a = a' mod qj. 



From Lemma|3.7l we have 



V- 1 loglog(g(V,g2) 



10l'i\a' 



and 



so 



(3.26) 



y 1 ^ log(gogi) ^ 
I « '70'? 1 



1 (loglog(q'o + ?i +^'2))^ 



iQi'l W 

q2\a—a' 



qhiq'm 



Substituting ( |3.26| l into ( |3.25| l gives the desired bound 



(3.27) 



□ 



Note that T] — 77 > since < T] < 1 .We may take 77 small enough so that ( 3. 12 1 and ( 3. 10 1 imply 



> (t7-ctj2)X>X 



as desired. 



□ 



A PROOF OF THE POSITIVE DENSITY CONJECTURE FOR INTEGER APOLLONIAN CIRCLE PACKINGS 19 

We note that the methods used here are easily generahzable to many discrete linear algebraic groups 
acting on with an integral orbit. If the group contains several Fuchsian subgroups as in the case 
of the Apollonian group, we may restrict to the orbits of these subgroups as in Section |2] We would 
again utilize the subgroup's preimage in the spin double cover of SO to relate the problem to integers 
represented by a binary quadratic form. This would yield a comparable lower bound on the number 
of integers less than X in the orbit of the group (counted without multiplicity). 
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